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چکیده بل‌ها یکی از مهم‌ترین سازه‌های هیدرولیکی ارتباطی به شمار می‌روند. برخورد جریان به پایةٌ پل‌ها فر سایش اطراف پایه ‏ وآ سیب‌دیدگی 
آن را به دنبال دارد. ازاین‌رو یکی از مباحث مهم در زمینة ایمنی پل‌هاء کنترل فرسایش بستر اطراف پایه‌ها است. تاکنون مدل‌های فیزیکی محتلفی در 
آزمایشگاه برای بررسی رفتار جریان و فرسایش بستر پیرامون آن سانحته شده‌اند که صرف‌نظر از هزینه و زمان‌بر بودن آزمایش‌هاء مشکلات ناشی از 
تغییرات مقیاس را نیز به همراه دا شته‌اند. امروزه با | ستفاده از کدهای دینامیک سییلات محا سباتی کا رآمدتر رفتار هیدرودینامیکی سیالات به‌طور 
گسترده‌تری بررسی می‌شود. در این تحقیق برای این منظور از نرم‌افزار ۳10۷310استفاده شده است. این پژوهش به شبیه سازی جریان اطراف مدل 
پایه‌های پل محتاف پردانعته و الگوی جریان و فرسایش بستر را برر سی می‌نماید. با توجه به حساسیت مل‌های عددی به تعداد سلول‌های شبکه 
در میدان محاسباتی» صخت سنجی تعداد محتلف شبکه با مقایسه نیمرخ‌های سرعت انجام یافته است. سه مدل آشفتگی مختلف نیز مورد بررسی 
قرا رگرفته و پس از برر سی نیمرخ‌های سط حآزاد آب الکوی جریان اطراف پایه و عدد استروهال, شبیه سازی جریان با مدل آ شفتگی ٩تلآانجام‏ 
کردیده است. تعداد ۱۵ مدل نیز برای ارزیابی دو پارامتر مهم عام لآب‌شسستگی» جدایش جریان و جریان رو به پایین, انتحاب شده‌اند. ابزارهای 
انتخابی در اين را ستاء سکی مقطع هند سی پایه و شکاف می‌با شند. نتایج دشان داد در اغلب موارد استفاده از یک ابزار کنترل آب شستگی بهتر از 
ت رکیب ابزارهای مختلف عمل می‌کند با به‌کارگیری شکل عدسی به‌عنوان مقطع هندسی پایه, عم قآب‌شستگی 7 9۸ درصد کاهش و با مدل ت رکیبی 
مقطع م ستطیل یگ ردگ و شه, سکو و شکاف بیشین عمق فر سایش یافته تا ۵۰ در صد کاهش ذشان داد. نتایج به‌د س تآمده بیانک رآن است که کنترل 
جدایش جریان و پدیدة شکست کردابه به‌مرائب بیش ا زکنترل جریان رو به پایین د رکاهش فرسایش بستر موثر است. 

واژه‌های کلیدی شبیه‌سازی عددی. رهایی متناوب گردابه‌هه نرم‌افزار 0 مدل آشفتگی. آب‌شستگی. پاية پل. 
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(۲) استاد. گروه مهندسی عمران. دانشکده مهندسی عمران دانشگاه صنعتی سهند. تبریز. 
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مقدامه 

اب‌ششتکی بلبلهاین استت که در آثر فرسایس پسی توسط 
جریان آب و حمل مواد بستر توسط نیرویی که این جریان 
به مواد بستر وارد می‌کند. به وجود می‌آید [1]. رادکیوی 
فا سکن را نی ی مرگ رود 
درنتیجه عمل فرسایش بستر توسط جریان و حمل مواد 
سخلا شمان ان توسط این جریان رخ می‌دهد [2]. سومر 
و فردسو بین آب‌شستگی و فرسایش کلّی. تفاوت قائل 
شده و آب‌شستگی را فرسایش مصالح بستر به سبب 
حضور یک سازه در معرض جریان دانسته‌اند [3]. 

مطالعات آزمایشگاهی فراوانی در زمينةٌ مدل‌سازی 
جریان و آب‌شستگی اطراف پایه‌های پل توسط 
پژوهشگران مختلف صورت گرفته است. به‌عنوان نمونه 
می‌توان به مطالعات لارسن و تاج شن و همکاران, ملویل 
و درگاهی. بروسرز و رادکیوی. گراف و ایسراتو 
دریسدل. گریمالدی و همکاران. گادیو و همکاران و 
گیه‌مو و همکاران اشاره نمود [4-12]. اغلب این مطالعات 
آزمایشگاهی مبتتی بر شرایط ایده‌ال‌سازی شده هستند 
لیکن اعمال شرایط پیچیده در مدل‌سازی آزمایشگاهی 
مشکل خواهد بود. ازاین‌رو یک فرایند رضایت‌بخش 
جهت بررسی جریان و انتقال رسوب. بهره‌گیری از مدل- 
های عددی است. نرم‌افزار 0 یک مدل عددی برای 
شبیه‌سازی سه‌بعدی جریان در هندسه‌های پیچیده است 
که توسط شرکت ۳۱0۹۵6866 در سال ۲۰۰۸ ارائه 
گردیده است. 

ریچاردسون و پانچانگ حفرة آب‌شستهُ ایجادشده 
در اطراف پاية پل استوانه‌ای را با ۳1030 مورد مطالعه 
قرار داده و نتایج غلادقن. تفس امه را با نتایج 
آزمایشگاهی ملویل و رادکیوی مقایسه نمودند. ایشان به 
اهمّیّت استفاده از اندازهٌ سلول‌های مناسب در شبکه‌بندی 
مدل تأکید کردند [13]. ین و همکاران جریان رو به پایین 
در بالادست پایه را در شرایط بستر صاف محاسبه نمودند 
تا با مقیاس‌بندی سرعت جریان در بالادست پایه شلّت 


جت جریان در حفرةٌ آب‌شسته را به دست آورند [14]. 


نشريه مهندسی عمران فردوسی 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


واسکو و والش به بررسی فرسایش بستر پیرامون پایه‌های 
مرکب با استفاده از اين نرم‌افزار پرداعتند. ایشان نشان 
دادند نتایج به‌دست‌آمده از, شییه‌سازی غددی اختلاف 
بسیار اندکی با نتایج تجربی موجود دارد [15]. عبدالعزیز 
و همکاران به مطالعةٌ عددی الگوی جریان و فرسایش 
بستر پیرامون پاية پل استوانه‌ای با نرم‌افزار ۲10۷39 
پرداختند [16]. مهرزاد به بررسی پدیدة آب‌شستگی حول 
پایه‌های مخروطی شکاف‌دار و بدون شکاف پرداخته و 
نشان داد که نتایج عددی و تجربی تطابق خوبی داشتند 
الگوی جریان در کانال مستقیم پرداخته و دریافتند انطباق 
نتایج بویت امه از شه‌سازی ‏ اغدهوی و نتایج 
آزمایشگاهی موید دقت بالای نرم‌افزار ۳۱0۷39 در شبیه - 
سازی الگوی جریان و گردابه‌های اطراف پایه است [۱۸]. 
حمیدی‌فر و امید (۱۳۹۱) شبیه‌سازی جریان سه‌بعدی در 
آب‌راهه‌های با مقطع مر کب را انجام دادند [۱۹]. توحیدی 
(۱۳۹۲) به‌صورت عددی به بررسی تأثیر پایه‌های با 
منحنی لگاریتمی در کاهش فرسایش بستر پرداخت. 
ایشان کاهش ۶۰ درصدی در عمق آب‌شستگی را گزارش 
نمود [۲۰]. حسن‌زاده و همکاران. تأثیر استفاده از عدل 
ترکیبی مقطع هندسی و شکاف را به‌صورت تجربی 
بررسی و نشان دادند به دلیل احتمال گرفتگی شکاف در 
زمان وقوع سیلاب. این ابزار نمی‌تواند گزينة مناسبی برای 
کت لش کباش ۲۲۱7 

سکو و شکاف میانی به‌صورت منفرد و نیز به‌صورت 
ترکیی جهت کترل. آب‌شستگی بيشنهاد. شده‌اند. .هداف 
از انتخاب این ابزارها مطالعه تأثیر دو پارامتر اصلی 
آب‌شستگی. جدایش جریان از وجوه جانبی و جریان رو 
به پایین. بر آغاز و گسترش آب‌شستگی است. در ترکیب 
سکو و شکاف. جریان رو به پایین کنترل شده و با 
به‌کارگیری مقطع هندسی پایه جدایش جریان از وجوه 
جانبی پایه کمینه می‌گردد. ترکیب این روش‌ها نیز در 
راستای کنترل توآمان دو پارامتر فوق بوده است. هر کدام 


سال سی و سوم شمارة دی ۱۳۹۹ 


نازیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


از ترکیب‌ها تأثیراتی در کنترل آب‌شستگی داشته و حتی 
برخی توانسته‌اند کاهش قابل توجهی در عمق 
آب‌شستگی ایجاد نمایند. بااین‌وجود دارای ضعف‌هایی 
بوده‌اند که شاید به لحاظ سازه‌ای و پا عوامل اجرایی برای 
زمان حال توجیه‌پذیر نباشند. این ترکیبات. شامل این 
مواردند: 
۱. مقاطع هندسی پایه (مقاطع دایره‌ای» مستطیلی 
گردگوشه و عدسی)؛ 
۲. پاية استوانه‌ای با سکو (ارتفاع‌های مختلف برای 
سکو بررسی شده است)؛ 
۳ ترکیب سکو مقطع هندسی پایه و شکاف میانی. 
بنابراین درمجموع تعداد ۱۵ مدل با نرم‌افزار 
0 مورد ارزیابی قرار گرفته‌اند. عمده‌ترین هدف 


تحقیق» بررسی دو عامل ایجادکنندة آب‌شستگی است. 


مواد و روش‌ها 
نم‌افزار 17101739 
جهت شبیه‌سازی عددی معادلات حاکم بر جریان‌های 
آشفته, از بستهٌ نرم‌افزاری ۳0۱۷39 استفاده شده است که 
توانایی محاسبات بالا در جریان کانال‌های باز را دارد. در 
این مدل. معادلات حاکم بر جریان با استفاده از تقریب- 
های حجم محدود حل می‌شوند. در این روش میدان 
محاسباتی به شبکه‌هایی با سلول‌های مستطیلی تقسیم- 
بندی شده که برای هر سلول میانگین کمیّت‌های وابسته 
وجود دارد؛ یعنی تمامی متغیرها در مرکز سلول محاسبه 
می‌شوند به‌جز سرعت که در مرکز وجوه سلول حساب 
می‌شود. در اين نرم‌افزار از دو تکنیک عددی حجم سیّال 
(۷0۴) و روش کسر مساحت- حجم مانع (۳۸۷۵8) 
برای شبیه‌سازی هندسی مرزها استفاده می‌شود. روش 
حجم سیّال برای نشان دادن رفتار سیّال در سطح آزاد آب 
(مرز آب‌وهوا) و روش کسر مساحت - حجم مانع برای 
شبیه‌سازی سطوح و مرزهای هندسی احجام صلب (مرز 
آب و جامد) کاربرد دارد. برای شبیه‌سازی اثرات 


آشفتگی. پنج مدل آشفتگی طول اختلاط پرانتل» مدل 


سال سی و سوم, شمارة دو» ۱۳۹۹ 


تک‌معادله‌ای. مدل دومعادله‌ای (2- 1 مدل 0( و 


مدل شبیه‌سازی گردابه‌های بزرگ در این نرم‌افزار معرفی 


شده است [22]. 


معادلات مورد استفاده در نرم‌افزار راز ۱ 
مدل هیدرودینامیک بر مبنای حل معادلات ناویراستوکس 
سه‌بعدی و معادلهٌ پیوستگی است. معادلات ناویراستوکس 
و شقات کین مورد استفاده در 71007310 برای جریان‌های 
غیرقابل تراکم در ادامه آمده است [22]: 


0 

)۱ 0 تس وج 
ی ی ی ادا بش 9۳1 
 ۲(‏ لا[ ی [ رنب + 


کت هنت شرمع ۳ فان کب شا یاه 
برای جریان در راستای :۷ کسر حجم باز برای جریان, 
تانب قاس از تتروهای بوشیی. 1 شا کاشین از 
نیروهای ویسکوزه1 ترم دراگ بین ذرات رسوب و م 
تسه سیال: استت» نش برتتی بپیعربا اتفاده: از تنم 
دیواره برای جریان‌های آشفته سه‌بعدی محاسبه می‌گردد 
[22,23]: 


۳( () 9 دا < ۱1 


که ا سرعت جریان. جل سرعت برشی» ۲ تنش 
برشی بستر» 6 فاصله از دیواره. ۷ وسیکوزیته 
سینماتیکی» ۶-۰۴ ضریب فن-کارمن و معا ضریب 
زبری نیکورادزه بت 


مدل آب‌شستگی 10۷380 مدل آب‌شستگی نرم‌افزار 
0 با استفاده از تقریب حجمی بقای جرم و معادلهة 
اتفالد ریا بش اسان مر ال فان افتال 
رسوب و آب‌شستگی بستر را برآورد می‌کند و با استفاده 
از روش کسر حجم سیّال ارتفاع سطح بستر در هر سلول 
مات زا پش‌تتی ی نماید: مدل اس‌شستکی ای 
نرم‌افزار شامل انتقال رسوبات بستر و رسوبات معلق 


نشریه مهندسی عمران فردوسی 


است. رسوبات بستر که توسط ذرات محصور شده‌اند 
به‌راحتی نمی‌توانند جابه‌جا شوند و درصورتی‌که 
فرسایش یابند و در سطح مشترک سیّال و رسوب 
ب‌صورت معلّق درآیند می‌توانند حرکت کنند. با فرض 
جریان دو فازی ذرات رسوب در سیّال و برقراری رابطهة 
موستم در فاز پپوسه سیال» سعاخلة (6) ورفاز تاپیوسته 
دای ترسوسی الق رها از تقایل شاد تفر )ان ماه 
(۵)» رابطة (3) حاصل می‌شود [22]. 


(( ۲ + 10۳- ۳0 + 2 
0 8 
گر _ تست ونا0 
(6 لا 1 ۲ + ۷۲ ولا ۰۷ + 
نگل - 1 01 - ر عزیونا0 

() بدا - نز + (7 9 ( و۳ . + انله۳۳ 


که‌یلا سرعت ذرات رسوب. + «(:1-1 <) تا 
(ونای سرعت متوسط حجمی سیال. نیروهای 
وزنی»(۷ - ولا <) لا سرعت نسبی ذرات رسوب و 
سیّال (ا - ولا <) »ولا سرعت رانشی ذرات رسوب 
است. رابطة بقای جرم برای جریان با سرعت متوسط 
حجمی سیال ۷۸-0 نیز باید برقرار باشد. تتش برشی 
سبب فرسایش و جابه‌جا شدن رسوبات در سطح بستر 
می‌گردد. این فرسایش تابعی از تتش برشی بحرانی ,7 و 
دانسیته سیّال و رسوب است. پارامتر شیلدز بحرانی .0 
کمینه تتش برشی جهت بلند کردن ذرات رسوب از سطح 
مرک الق نس هقال را نشان نهد 124 


۷ رتم۲ ۳ 0 
5 سوم ۳ هه 


مه قطر ذره رسوب و اندیس 2 به ۲ امین ذره 
رسوب اشاره دارد. مقدار پایه برای عدد شیلدز بحرانی 
در بسترهای صاف افقی تعریف و معمولاً برابر مقدار 
ثابت ۰/۰۵ در نظر گرفته شده و یا از رابطةٌ سالبوی و 
وایت‌هاوس مقدار دقیق آن تعریف می‌گردد [25]: 


۸0( رد0020 - 0.0551 + سش. 2 ورب9 


وب 1+1.2 


5۳ اندازه بی‌بعد 1 امین ذره رسوت است. 
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شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


در این مدل پدیدهٌ فرسایش و رسوب‌گذاری هر دو 
پدیده‌های متقابل میکرو هستند که به‌صورت هم‌زمان رخ 
می‌دهند. ترکیب این دو پدیده به‌منظور بررسی تغییرات 
حالص بین رسوب بستر و رسوب معلّق انجام می‌گيرد. 
برای محاسبهٌ فرسایش بست سرعت بلند شدن و کنده 
شدن هر ذرّه رسوب از سطح بستر» «روزلا که از سرعت 
برشی بحرانی بستر بیشتر است بر اساس معادلة وینرپ 
و همکاران (۱۹۹۲) حاصل می‌شود: 


(1 - وه) م20 روش 9 08 3 جر و101 < صعنبلا 
(٩)‏ 
مه ضریب بلند شدن هر ذرّه رسوب و و بردار نرمال بر 
سطح بستر رسوبات است. در محاسبة ته‌نشینی رسوبات. 
سرعت ته‌نشينر ارائه شده توسط سالبوی (۱۹۹۷) مورد 
استفاده قرار ین که: 
۷۶ 1/2 1 
|1036 - ( 104903 + 0 < :۱0:2 
(۱۰ 
که 2 شتاب تقل و طهلتهعلا هم‌راستا با شتاب ثقل 
برای ذرّه رسوب « به‌صورت زیر تعریف می‌شود: 
)۱۱ 1 هرا < و 
1[7-:«2)5] 
م0 نرخ انتقال رسوب حجمی به ازای واحد 
عرض بستر و « بر اساس معادله مییرپیتر-مولر (۱۹۸۶) 
محاسبه می‌شود [26]: 
)۱۲( (ت 9 - و«0)وظ < و9 


8 ضریب بار بستر بوده و مقدار آن برای نرخ 
انتقال رسوب پایین بین ۵ تا ۵/۷ برای نرخ انتقال متوسط 
حدود ۸ و برای نرخ انتقال بالا برابر ۱۳ اختیار می‌گردد. 
«و6 جزء حجمی ذرّه رسوب « در بستر است. رابطةً 
ون-رایان (۱۹۸۶) برای محاسبهً ضخامت لایة رسوبی م9 


مورد استفاده قرار می گیرد [27]: 


سال سی و سوم شمارة دی ۱۳۹۹ 


نا زیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


۳ 15 - س) و1 0.30 ع ورط 
هندسة مدل 
ملویل (۱۹۷۵) استفاده گردیده و کانال مورد آزمایش در 
پژوهش ملویل با طول ۱ متر. عرض ۰/۶۵۱ متر و عمق 
جریان ۰/۱۵ متر شبیه‌سازی شده است. دبی جریان 
۲ لیتر بر ثانیه میانگین سرعت جریان ۰/۲۵ متر بر 
ثانیه» قطر پایه استوانه‌ای ۵/۱ سانتی‌مت متوسط قطر 
رسوبات بستر ۰/۱۳۸۵ میلی متر؛ زاویه ایستائی ذرات 
رسوب ۳۲ درجه و چگال نسبی آن‌ها ۲/۷۸ است [6]. در 
مطالعهٌ عددی پاية دایره‌ای در وسط کانال تعبیه شده است 
تا در این فاصله جریان قبل از رسیدن به پایه به حالت 


توسعه‌یافته رسیده باشد. 


صحت سنجحی نتایج عددی 
تعریف شرایط مرزی 
در بلوک محاسباتی به‌طور جداگانه شرایط مرزی شش- 
گانه تعریف می‌گردد. برای دیواره‌ها شرط مرزی غیر 
لغزشی دیواره» برای ورودی کانال شرط مرزی فشار معیّن 
با اعمال تراز سطح آب و سرعت جریان در راستای " 
شرط مرزی جریان خروجی برای خروجی کانال و برای 
سطح آزاد جریان شرط مرزی تقارن در نظر گرفته شده 


است. 


بررسی انداز؟ سلول‌های میدان محاسباتی 
حساسیت‌های مدل عددی به شبکه‌بندی و انفصال میدان. 
همواره یکی از مسائل مهم در مدل‌سازی عددی بوده 
است. برای مشخص کردن حدود شبکه‌بندی. بلوک‌های 
مستطیلی شکل در نرم‌افزار ۳30۷3 ارائه شده است که 
کليهٌ اندازه‌های سازهٌ مورد نظر و فضای اشغال‌شده توسط 


سیّالات در داخل آن‌ها تعریف می‌شود. برای مدل عددی 


سال سی و سوم. شمارهُ دو, ۱۳۹۹ 


حاضر از یک بلوک شبکه‌بندی استفاده شده و تعداد 
سلول‌های محاسباتی مختلف مورد آزمون واقع شده‌اند. 
در محدودهٌ اطراف پایه نیز به دلیل اینکه تمامی 
برجستگی‌های سازه توسط مدل عددی شناخته شوند و 
هیدرودینامیک مدل با دقّت بیشتری شبیه‌سازی گردد نیاز 
به استفاده از شبکه‌بندی ریزتر است. تعداد سلول‌ها و 
استفاده از سلول‌های ریز در اطراف پایه در شکل (۱) 


نشان داده شده است. 


شکل (۱): تعداد سلول‌ها در نقاط مختلف شبکه‌بندی میدان 


محاسباتی 


به‌منظور انتخاب تعداد سلول مناسب. مدل‌سازی 
تغییرات سطح آزاد آب در وسط کانال با استفاده از چهار 
تعداد مختلف سلول. با مقادیر ۱۳۰۰۰۰ ۲۵۰۰۰۰ 
۰ و 10۰۰۰۰ انجام يافته و نتایج به‌دست‌آمده با 
نتایج تجربی ملویل (۱۹۷۵) مقایسه شد (شکل ۲). در 
شبیه‌سازی حاضر مدل آشفتگی 15 به‌کار برده شده 
است. در جدول (۱) ارتفاع محاسبه‌شده سطح آب توسط 
نرم‌افزار با داده‌های تجربی ملویل (۱۹۷۵) مقایسه شده و 
برای هرکدام از شبکه‌بندی‌هاء مقدار خطا در مقطع عرضی 
معیّن شده محاسبه شده است. برای محاسبةٌ خطا از رابطة 


زیر استفاده شده است: 
ره 0 4 ۳ ۴ - 8 


که در آن ز۷ عمق تجربی سطح آب ور عمق 


جریان تفس آسله از ییاز لو استه: 


نشریه مهندسی عمران فردوسی 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


وت 


(1975) ۲۷۲۵۱۷۱1۱6 -امصمصدزه مور ۰ 


۷/65 130000 اممزرموریا]( سس 
۷6۶ 250000 -1هعزع هب۲ 
5 360000 م۲۵۲6 


تراز سطح آب (70) 


طول کانال (7) 


شکل (۲): مقایسةُ نیمرخ طولی سطح آب محاسباتی با داده‌های تجربی ملویل (۱۹۷۵) برای تعداد سلول: 
الف) ۰ ب) چا ۱۲۳۹ 


جدول (۱): مقايسة ارتفاع محاسباتی سطح آب با داده‌های تجربی ملویل (۱۹۷۵) برای شبکه‌های مختلف 


تعداد شبکه 
شبکه‌بندی درشت 
شبکه‌بندی متوسط 
شبکه‌بندی ریز 


شبکه‌بندی بسیار ریز 


از نتایجم جدول فوق می‌توان دریافت خطای 
شبکه‌بندی درشت نسبت به دو شبکه‌بندی متوسط و ریز 
بسیار زیاد است که درنتیجه برای سایر شبیه‌سازی‌ها؛ از 
آن صرف‌نظر می‌شود. از طرفی میزان بهبود خطا در 
شبکه‌بندی ریز نسبت به شبکه‌بندی متوسط بسیار 
این موضوع نشان می‌دهد با افزایش تعداد شبکه به 
مقداری بیش از آن و افزایش حجم سلول‌های محاسباتی, 
مقدار جواب مستقل از تعداد سلول‌ها باقی می‌ماند. پس 
تعداد ۳۱۳۰۰۰ سلول با اطمینان غویی برای: سار مدل‌ها 


مورد استفاده قرار گرفته است. 
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خحطای 5 () 
۱۳۹۰ ۱-۱( 
۱۳۵۰۰۰ 9۱۹-۶ 
۳ ۱/۳۹۰۸ 
...1-9 ۱/۸۳۳ 


بررسی مدل‌های آشفتگی 

منحنی‌های سرعت جریان 
پس از انتخاب شبکه‌بندی مناسب. بررسی منحنی سرعت 
و الگوی جریان برای انتخاب مدل آشفتگی انجام شده 
است. در شبیه‌سازی عددی جریان‌های ثئانویَةُ پیرامون 
پایه» اعمال اثر آشفتگی بسیار حائز اهمَیّت است؛ ازاین‌رو 
سه مدل آشفتگی 2 - 1 استاندارده 116 و 15 انتخاب 
گردیده و مورد آزمون قرار گرفته‌اند. شکل (۳) مقايسة 
نتایج عددی و تجربی منحنی‌های سرعت در صفحه 17 
را نمایش می‌دهد. این منحنی‌ها در فاصلة د۵ در 
بالادست محل پایه در نظر گرفته شده‌اند که 0 قطر پایه 


اشتتت: 


سال سی و سوم شمارة دی ۱۳۹۹ 


نا زیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


کااعع؟ ام)معصننهیرظ وعااز۷ ۷۵۱( 


۵۵۱ هوعاباطا:۲ ۲۵۱۵ - لهءزتعصصب( 
۱۷۵0۵1 عمصعلناطاین ۲ فظ - لممزرمریا]( سس 


۱۹ 


سرعت در راستای * (7/5) 
شکل (۲): مقايسة منحنی‌های سرعت حاصل‌شده با سه مدل 


آشفتگی با منحنی سرعت مدل تجربی ملویل (۱۹۷۵) 


برای منحنی به‌دست‌آمده از مدل آشفتگیه - 1 
استاندارد مقدار خطا برابر ۸/۱۰۹ درصد با مدل آشفتگی 
0 برابر 1/۵76 درصد و با استفاده از مدل کلب پرابر 
۲ است. مشاهده می‌شود استفاده از مدل آشفتگی 
5 موجب نتایج دقیق‌تر شده و ازاین‌رو در سایر 
شبیه‌سازی‌ها مورد استفاده قرار گرفته است. شکل (4) 
مقایسهة سرعت جریان محاسباتی با داده‌های تجربی 
ملویل را در نقاط مختلف بستر نمایش می‌دهد. 

تطابق قابل قبول مقادیر محاسباتی و تجربی سرعت 
برای مدل آشفتگی 1,55 تا اندازة زیادی در شکل مشخص 
است. شبیه‌سازی دقیق جریان به‌ویژه در نزدیکی مرز پایه 
بسیار حائز اهمَیّت است زیرا نقش قابل توجهی در انتقال 
بار بستر و فرسایش بستر پیرامون پایه خواهد داشت. 


03 


۸ 1:2 021 
۲۲ 13۵6 1 
6 2 


سرعت اندازه گیری شده (۱7/5) 
ت‌‌ 


05 0 05 
سرعت پیش بینی شده (۷5) 


شکل (۶): مقايسه مقادیر محاسباتی و تجربی سرعت در نقاط 
مختلف پیرامون پایه با استفاده از سه مدل آشفتگی 


سال سی و سوم. شمارهُ دو, ۱۳۹۹ 


بررسی الگوی جریان 

برای ارزیابی بیشتر مدل‌های آشفتگی در نرم‌افزار 
0 کانتورهای سرعت جریان و گردابه‌های 
برخاستگی ناحية دنبالةٍ پایه مورد مقایسه قرار گرفته‌اند. 
شکل (۵) کانتورهای به‌دست‌آمده با مدل آشفتگی 8 - 1 
را در تراز بستر نمایش می‌دهد. برای بررسی تناوب 
گرذانه‌شا: اتعرای مدل 1۲ *۸ کانبه به طول اتجامیده اسق: 

در شکل (1) گردابه به‌دست‌آمده از شبیه‌سازی با 
مدل آشفتگی ع - نشان داده شده است. در اعداد 
رینولدز بسیار کوچک ( 2 >18) جریان اطراف پایه‌ها کاملا 
به پایه چسبیده و جدایش جریان از مرزهای پایه رخ 


نمی‌دهد. 


شکل (۵): کنتورهای سرعت جریان در راستای طولی با مدل 
آشفتگی» - 6 در بازه‌های زمانی ۷ ثانیه‌ای 

با افزايش عدد رینولدز. جریان از لبه‌های پایة جدا 
شده و یک جفت گردابه متقارن در پشت پایه به وجود 
می‌آید. با افزايش عدد رینولدز, گردابه‌های متقارن پشت 
پایه از حالت دائم خارج شده و به‌صورت نوسانی درآمده 
و به‌طور متناوب در جریان پخش می‌شوند. در این مدل. 
عدد رینولدز برای جریان عبوری از اطراف پایه. برابر 
2-۲۳ است که با توجّه به بالا بودن عدد رینولدز, 
تشکیل گردابه‌های متقارن و عدم پخش آن‌ها در جریان 
منطقی نیست. مشاهده می‌ گردد که در زمان‌های ۰ و ۷۲ 


نشریهٌ مهندسی عمران فردوسی 


ثانیه. گردابه‌های برخاستگی کاملاً متقارن باقی‌مانده و 
نوسان گردابه در آن‌ها رخ نمی‌دهد. 

حالت متقارن گردابه‌ها و عدم پخش آن‌ها در 
جریان نشان می‌دهد این مدل آشفتگی نمی‌تواند شکست 
گردابه‌ها در پشت پایه را به‌درستی شبیه‌سازی نماید. 
شکل (۷) کانتورهای سرعت جریان در تراز بستر و 
گردابه‌های جریان شبیه‌سازی‌شده با مدل آشفتگی 110 
را نمایش می‌دهد. 

در این مدل نیز تغییر قابل توجهی در کانتورهای 
سرعت نسبت به مدل آشفتگی 2 - 16 حاصل نشده است» 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


لیکن نوسان گردابه‌ها در پشت پایه. در مقایسه با شکل 
() محسوس است. در مدل‌سازی عددی با مدل آشفتگی 
6 گردابه‌ها در زمان ۱۰ و ۲۰ ثانیه حالت متقارن 
دارند. در زمان 1۷ ثانیه» حالت متقارن گردابه‌ها حذف 
گردیده و رشد و گسترش گردابه‌ها در کناره‌های پایه 
مشاهده می‌ گردد. ابعاد گردابه‌ها در زمان ۷۶ انیه نست 
به زمان 1۷ ثانیه کاهش یافته است لیکن در وجوه جانبی 
پایه نوسانی مشاهده نمی‌گردد. همچنین فرآیند انتشار 
گردابه‌ها نیز درست نبوده و ناحيةٌ برخاستگی گردابه‌ها 
تنها به فاصلهٌ کمی از پایین‌دست پایه محدود شده است. 


شکل (۱): فرم متقارن گردابه برخاستگی در ناحیهٌ دنبالٌ جریان شبیه‌سازی‌شده با مدل آشفتگی ع - 16 


۲۱۵۱۲ ۱۵۲60۵ 


"سس 


و سح 
1 


۰ 
ِ 
۳ 


شکل (۷: (الف) کانتورهای سرعت جریان شبیه‌سازی‌شده و (ب) گردابه برخاستگی در ناحیهٌ دنبالٌ جریان با مدل آشفتگی 60 در 
زمان‌های مختلف 


شکل (۸): کانتورهای سرعت جریان شبیه‌سازی‌شده با مدل آشفتگی ور 


نشریه مهندسی عمران فردوسی 


سال سی و سوم شماره دی ۱۳۹۹ 


نازیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


در شکل (۸) کانتورهای سرعت جریان در تراز 
بستر و یک سیکل شماتیک از مراحل ایجاد. رشد و پخش 


نشان داده شده است. 

شکل بالا نشان می‌دهد گردابه‌ها در مدل آشفته 
8 به‌صورت نوسانی شبیه‌سازی شده و وقوع پدیده 
رهایی متناوب گردابه‌ها مشاهده می‌گردد. در فاز ۱ 
نوسان. گردابه کوچکی در وجه بالایی پایه ایجاد شده و 
گردابه بزرگ وجه پایینی از مرز پایه جدا شده و در حال 
پخش در جریان است. در فاز ۲ و با گذشت ۲ ثانیه از 
فاز ۱. گردابه کوچک بالایی رشد يافته و در حال پخش 
به سمت پایین‌دست است. هم‌زمان با رشد و انتشار 
گردابه بالایی به پایین‌دست. گردابه کوچکی در وجه 
پایینی در حال رشد کردن بوده تا اينکه به بیشینه توان 
خود برسد. در فاز ۳ گردابه بزرگ بالایی کاملا تشن 
شده و گردابه پایینی رشد پیدا کرده است. در فاز ۶ 
کر هناش ی ان ون به تیگ پاسی ر نصا سار 
بوده و گردابه بالایی در حال تشکیل است. درنهایت در 
فاز ۵ و پس از گذشت ۷ ثانیه. الگوی یکسانی از رشد و 
گسترش گردابه‌ها (در فاز ۱) مشاهده می‌گردد. فرکانس 
نوسان این گردابه‌ها (8). در حدود ۰/۸۹۷ هرتز است. با 
شروع پدیده نوسان گردابه‌ه نیروهای نوسانی براً و پسا 
از طرف جریان بر پایه وارد می‌گردد. ماهیّت نوسانی این 
نیروها به‌گونه‌ای است که می‌توان عملکرد این نیروها را 
در یک سیکل کامل ۷ ثانیه‌ای (دوره تناوب 1) منطبق بر 
سیکل نوسان گردابه‌ها در نظر گرفت. 

یکی دیگر از مشخصه‌های جریان غیردائمی. عدد 
استروهال جریان است. در پدیده رهایی متناوب گردابه- 
هاء عدد استروهال به‌عنوان اندازه فرکانس بی‌بعد پدیدة 
شکست گردابه در نظر گرفته می‌شود. عدد استروهال به- 
صورت 80/۷-<5 تعریف می‌شود که 56 عدد 
ابعروهان. خر کاشی گرفانه‌ها:ظ ط. باه و 0۳ 
سرعت جریان است. عدد استروهال به عدد رینولدز 
جریان و تعداد شبکه بسیار حساس بوده و با تغییر در 


مقدار آن‌ها؛ دچار تغییر می‌گردد. طبق نتایج عبداللهی و 


سال سی و سوم شمارهة دو» ۱۳۹۹ 


عاطفی (۱۳۹۰), در محدوده ۳۵۰> ۰۱۵۰6 حساسیت 
نتایج محاسبات به نوع شبکه‌بندی وجود دارد اما با 
افزایش عدد رینولدز و در محدوده ۳۵۰< 66 عدد 
استروهال مستقل از تعداد سلول‌های شبکه باقی می‌ماند. 
علّت تأثیر قابل توجه نوع شبکه‌بندی بر نتایج مربوط به 
عدد استروهال در محدوده ۳۵۰> ۱۵۰۲۵ را می‌توان به 
تغییر الگوی جریان و نفوذ جریان‌های گردابه‌ای از پایین - 
دست پایه به سمت لبه‌های جلویی, مربوط دانست [۲۸]. 
در این پژوهش با توجه به بالا بودن عدد رینولدز جریان 
(۰۳ ۱ عدد استروهال مستقل از عدد رینولدز و تعداد 
سلول‌های شبکه است. مقايسة فرکانس نوسان جریان 
گردابه‌ای و عدد استروهال برای نتایج تجربی ملویل و 
نتایج عددی حاضر در جدول (۲) آمده است. 


جدول (۲): مقايسة مقادیر تجربی و عددی پارامترهای 


هیدرودینامیکی جریان 


پارامتر نماد تجربی عددی 

دوره تناوب (ثانیه) 1 ۹/۳۷ ۷ 
فرکانس (هرتز) 0 ۰/۹۸۵ ۰,۸۹۷ 
عدد استروهال ۹ ۰/۲ ۰۱۸۳۳ 


عمق آب‌شستگی بستر 
عمق آب‌شستگی یک ساعت اول با توجه به نتایج تجربی 
ملویل (۱۹۷9) ضکت‌سنجی شده.و سیس یون تغییر دز 
پارامترهای رسوب مدل عددی. سایر مدل‌سازی‌ها انجام 
یافته‌اند. پارامتر شیلدز بحرانی. ضریب دراگ ضریب 
انتقال ذرات. ضریب بار بستر و زاوية ایستایی از مهم‌ترین 
پارامترهای رسوب به شمار می‌روند که باید در تنظیمات 
الیٌ مقدار صحیح آن‌ها از طریق سعی و خطا مشخص 
گردد. محدوده مناسب برای عدد شیلدز بحرانی 
0/02-1 است که بر حسب نمودار شیلدز. مقدار آن 
پرابر ۰/۰۳۶ به‌دست‌آمده است. مقدار ضریب دراگ برابر 


۵ برای ذرات کروی است. انگلاند و هانسون در سال 


نشریه مهندسی عمران فردوسی 


۷ پیشنهاد کردند که برای شن و ماسه مقدار ۱/۵ 
مناسب است و ساير مقادیر معمول آزمایشگاهی بین 
۵ است. در این بررسی مقدار ضریب دراگ 
برابر ۱/۲ وارد شده است. ضریب بار بستر که در رابطهةً 
میرپیتر-مولر وجود دارد نرخ انتقال بار بستر به ازای تنش 
پزشی بیش از تن پرشین بحرانی است: معدار این متریپ 
پایین و مقدار ۱۳ برای شلّت بالای انتقال ماسه تعریف 
ضریب تعیین شده است. یکی از پارامترهای مهم در مدل 
آب‌شستگی -رسوب» زاوية ایستایی ذرات رسوت است. 
آزمایشگاهی آن یعنی ۲۲ درجه وارد شده است. 
در مطالعه آزمایشگاهی ملویل. مقدار عمق 
آب‌شستگی تعادلی برای مدل مبنا برابر ۱۵ سانتی‌متر و 
مق آب شست گرم ۱ ساعته معادل ۸ سانتی‌متر است. در 


جدول (۲) مقايسة داده‌های تجربی و نتایج عددی ارائه 


شده است. مقدار خطا در شبیه‌سازی یک‌ساعته براپر ۵/۶ 
درصد به‌دست‌آمده است. این موضوع نشان‌دهندة تطابق 
فان کول شاد کدی و معادیر هر سلیی انشت: 

برای بررسی و مقایسه کامل‌تر نتایج تجربی و 
عددی. از رابطة ارائه شده توسط ملویل و چیو (۱۹۹۹) 
استفاده شده است. رابطة تجربی ملویل و چیو (۱۹۹۹۵) 
به شکل زیر است [29] 


۱0۱ | |(۵) تا 


که زمان تعادل از رابطة زیر به‌دست می‌آید: 


رم ۱ / ظ 
6 < 309 
"0 2 0 [1 


9 04( 1 
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نشريه مهندسی عمران فردوسی 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


جدول (۲): مقايسة مقادیر تجربی و عددی عمق فرسایش در 


بازه‌های زمانی ۵ دقیقه‌ای 


عمق آب‌شستگی (صطه) 
زمان (دقیقه) 
عددی تجربی 
۱/٩۱ ۵‏ ۱ 
۳ ۱/۳۷ ۱/۵ 
۳ ۱۳:۷ ۱/۳۹۰ 
۳۰ ۱/۳۹۹ ۱/۳/۸۰ 
۲۵ ۳/۱ ش 
۳۰ ۳ ۳/۲ 
۳۵ ۳/۷۰ ۳/۹ 
۰ 2-۱۱ 1/۲ 
1۵ 0/05۸ 0/۸ 
0۰ ۰134 1/۹۰ 
0 ۷/۳۸ ۷/۶ 
5 ۷/۹۲ ۸ 
خطای متوسط کل (/) :1 


در شکل )٩(‏ مقایسه نتایج عددی حاصل از شبیه- 
سازی و نیز داده‌های به‌دست‌آمده از رابطهٌ (۱۵) نمایش 


داده شده است. 


شکل (4): مقايسهٌ داده‌های حاصل از رابطة ملویل و چیو (۱1۹۹۹) 


و نتایج عددی 


مدل یایه‌ها 
پس از اطمینان از صحّت نتایج هیدرودینامیکی و رسوبی. 
سایر مدل‌های انتخابی مورد بررسی قرار می‌گيرند. مدل- 
های پیشنهادی در چهار گروه کلّی طبقه‌بندی می‌شوند. 


سال سی و سوم شمارة دی ۱۳۹۹ 


نازیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


گروه اوّل (۸) پایه‌های با مقاطع هندسی مختلف (دایره- 
ای. مستطیلی گرد گوشه بیضی و عدسی) و گروه دوم 
(8) پایه‌های سکویی با چهار ارتفاع ۵ 7 ۷و ۸سانتی‌متر 
برای ارتفاع سکو را شامل می‌شوند. اين اعداد از بررسی 
پیاشهاد بعرهیگران مصلت اععات ده است ,در کروه 
سوم (گروه 6)؛ مدل مرکب سکو و مقطع هندسی پایه و 
در گروه (1 ترکیب مدل‌های گروه (6) با شکاف بررسی 
شده است. در تمامی مدل‌هاء پاية استوانه‌ای با مقطع دایره 
به‌عنوان مدل مبنا بوده و بیشینهةٌ فرسایش ۱ ساعته بستر 
در بالادست مدل‌ها. به دست آمده است. پایه‌های مورد 
بررسی دارای عرض ۶ سانتی‌متر نسبت طول به عرض 
۲ متوسط قطر رسوبات ۰/۷۷ میلی‌متر: سرعث جریان 
۵ متر بر انیه عمق جریان ۲۰ سانتی‌متر» زاوية 
ایستائی ذرات ۳۲ درجه و چکالی نسبی ذرات ۲/۱۷ 
است. 

گروه اوّل» شامل مقاطع هندسی دایره مستطیل گرد 
گوشه. بیضی و عدسی است (شکل ۱۰-لف). در گروه 
دوم پایه‌های سکویی دارای سکویی به قطر 1 سانتی‌متر 
بوده (از هر طرف ۱ سانتی‌متر بیرون آمدگی) و ارتفاع ۵ 
1 ۷ و ۸ سانتی‌متر نسبت به تراز بستر می‌باشند (شکل 
۰-ب). در پایه‌های شکاف‌دار نیز عرض شکاف ۲ 
سانتی‌متر و ارتفاع آن با توجّه به ارتفاع سکو در بالای 
تراز بستر تعیین شده است (شکل ۱۰-پ). در جدول 
(4) و شکل (۱۱) جزئیات مدل پایه‌ها ارائه گردیده است. 


۱۹۹۱ اتدامما1 1 یدای ۱۷ ۳ 

(الف) 
اج 4 
1 ایا اک ۱۲۵ 
۳ 
سب 

1 2۱ 
۱ ۱ ‌ 1 اوسا تما[ ۹ 
ات 
اما اص۲) 
(ب) 


سال سی و سوم شماره دی ۱۳۹۹ 


(پ) 


شکل (۱۱): (الف) مدل‌های گروه ۸ (ب) مدل‌های سکویی گروه 
3 (پ) مدل پایه‌های شکاف‌دار گروه 1 


1 
مثله مفله مسلشن 
سم ۳۹ 
تس 
--. مه ت سه 2 5 
سح(] ءس»(] 
مه هه ف‌ "۱ 
۱ ک وب سم "۱ 
۰ ۳1 0 
۲ ۳۵ ۱۱ ِ ۰ ۰ 4 
۳ 0 ۳7 
۳ سه امه 
‌ُ ‌ 8 ...اس ۳ 


شکل (۱۱): ابعاد مدل پایه‌ها 


در شبیه‌سازی این گروه از مدل‌هاء تعداد سلول‌های 
شبکه ۳۲۰۰۰۰ مدل آشفتگی 185 و زاوية ایستائی ذرات 


۲ درجه در نظر گرفته شده است. 


نتایج و بحث 
مدل میا 

در شکل (۱۲) فرسایش بستر و بردارهای سرعت جریان 
پیرامون مدل مبنا نشان داده شده است. 

مشاهده می‌شود جریان نزدیک شونده پس از 
برخورد به پایه به سمت بستر منحرف شده و جت جریان 
در مجاورت آن تشکیل شده است. با برخورد جت جریان 
به بستر بالادست. گردابه‌های نعلاسبی تشکیل شده و در 
پایین‌دست پایه نیز بازگشت جریان و تشکیل گردابه‌های 
برخاستگی مشاهده می‌گردد. 


نشریه مهندسی عمران فردوسی 


۳۲ شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


براق هل متا پیفینه عمق. آب‌شستگی. رک‌ساهته می/(0 به کار برده شده است؛ که.»4 بیشینة عمق 
برابر ۲۹ میلی‌متر است. برای محاسبهٌ درصد کاهش آب‌شستگی پایه مبنا و وه بيشينةٌ عمق آب‌شستگی 
آلش تس کزن مدل‌های پیشنهادی. رابطةً ب م0 و1۲0 مدل‌های مورد بررسی اسسیت: 


انحراف جریان رو به پادین 


تشکیل گردابه‌های نعل اسبی 


یت در بالادست پابه 


شکل (۱۲): فرسایش بستر بالادست و پایین‌دست پایه به همراه تشکیل گردابه‌های نعلاسبی و برخاستگی 


جدول (4): مشخصات مدل پایه‌های پل 


نام احتصاری 
گرره | شیاره مد | ۲ سق ثِِ« هم | هه اه اه و۵ 

مدل مین ۱ 0 دایره‌ای 1 7 ات | ]۳ 
۲ تن مستطیلی گرد گوشه ۴ | تم: | سس | ۲۳۳۵۶ 
۵ ۳ 3 بیضی 1 1 ست | ات | ات | و۳ 
3 1 عدسی 3 1 ت ت [ ۳ 
۱ 051 دایره‌ای با سکوی ۵ سانتی‌متری 3 3 1 ۵ تست | ۳ 
ِ 092 دایره‌ای با سکوی 1 سانتی‌متری 3 5 4 ٩‏ | جب. | معا 
۱ ِ 083 دایره‌ای با سکوی ۷ سانتی‌متری 3 1 1 ۷ بر ۰| ۳9 
‌" 00۹4 دایره‌ای با سکوی ۸ سانتی‌متری 3 4 5 | عت. |« 
1 3 مستطیلی گرد گوشه با سکو ۶ | ٩‏ | ۱ | ؟ | -- |۳۰ 
0 ۰ ۳5 بیضی با سکو 3 1 1 1 | ۳۹ 
۷۱ 113 عدسی با سکو 3 ۹ ۹ ۴ | بت | م۳ 
۲ 55 دایره‌ای با سکو و شکاف میانی 3 1 7 7 ۲ ۳۹۰۰ 
۳ 155 یی 7 1 1 1 ۲ ]| ۳۹۰۰ 

17 شکاف میانی 
۱ ۱ بیضی با سکو و شکاف میانی 3 1 1 1 ۲ ۳۹۰۰ 
۱۵ 11۳۹ عدسی با سکو و شکاف میانی 3 1 1 1 ۲ ۳۹۰۰ 


تقمر 4 متسین مرن فرقارستی سال سی و سوم شمارة دی ۱۳۹۹ 


نازیلا کاردال- حبیب حکیم زاده- یوسف حسن زاده 


مدل‌های گروه ۸ 

آشفتگی در آرایش خطوط جریان در پایه‌های با دماغة 
تیز و آیرودینامیک کمتر از حالت دماغةٌ پهن بوده و شکل 
پایه تأثیر قابل ملاحظه‌ای در حذف گردابه‌های نعل‌اسبی 
و برخاستگی پیرامون پایه دارد. در اين راستا سه شکل 
پایه مستطیلی گرد گوشه. بیضی و عدسی که به ترتیب 
میزان انطباق آن‌ها با الگوی جریان بیشتر شده بررسی 
شده‌اند. در جدول (0) نتایج حاصل از شبیه‌سازی 
یک‌سافته برای مثل‌های کر وه ۸ آرانه شنده اس 


جدول (۵): حداکثر عمق آب‌شستگی مدل پایه‌های گروه ۸ 


مدل 9 ۱9/۳ 4 13 
۳/0 
پایه | (صه) | ره | (صص) 2 
6 3 3 ۲۹ ۵ اس 
8 3 3 ۲۱ ۰/۳۷۵ ۲۷ 
3 ۶۰ ۶۰ ۱3 ۰/۱۳۳۵ 31 
۳ ۰ 3 ۱۲ 2 5۸ 


در شکل (۱۳) و (۱۶) نمایی دوبعدی و سه‌بعدی 
از فرسایش بستر پایه‌های با مقطع مستطیلی گردگوشه و 
عدسی نشان داده شده آنتخنت همان‌طور که از شکل‌ها 
معلوم است فرسایش بستر اطراف مدل با بسیار کمتر از 
مدل ٩‏ است. 

با استفاده از مقطع هندسی عدسی. آب‌شستگی تا 
19۸ در مقطع بیضی تا 16 و با مقطع مستطیل گرد 
گوشه. آب‌شستگی تا ۸۲۷ کاهش یافته است. این کاهش 
نشان می‌دهد مقطع عدسی بهتر از سایر مقاطع عمل کرده 
و تأثیر بیشتری در کنترل فرسایش بستر داشته است. تأثیر 
مقطع بیضی نیز در کاهش آب‌شستگی قابل ملاحظه 
است. در مقطع عدسی به دلیل انطباق زیاد مقطع با الگوی 
جریان. الگوی جریان کمتر دستخوش تغییر شده و 
جداشد گی جریان تنها محدود به یک نقطه می‌گردد؛ 
بنابراین جداشد گی جریان از پایه حذف گردیده و احتمال 
تشکیل گردابه‌های نعل‌اسبی نیز در بالادست پایه ضعیف 
می‌گردد. 


درنتیجه با تضعیف پارامترهای عامل 


سال سی و سوم شماره دی ۱۳۹۹ 


آب‌شستگی (جدایش جریان و جت جریان رو به پایین) 
به میزان قابل توجهی از شدّت آب‌شستگی کاسته می- 
گرد 


09 
با 


اه ] اقا قاتا د. 


9 
1 
۳ 


(الف) (ب) 


(الف) (ب) 


شکل (۱۵) تغییرات آپ‌شستگی یک‌ساعته بیشینه 
را برای مدل‌های گروه ۸ نشان می‌دهد. در مدل‌های ) و 
8 تغییرات فرسایش بستر در انتهای یک ساعت قابل 
ملاحظه است؛ لیکن در مدل‌های 8 و ب تمایل برای 
رسیدن به تعادل در زمان یک ساعت مشاهده می‌شود. به 
عبارتی نمودار فرسایش در این مدل‌ها سریع‌تر به مجانب 
فقی میل کرده و تغیبرات فرسایش بستر در این مدل‌ها 
تنها در لحظات اولیه قابل ملاحظه است. 


مدل‌های گروه 8 

در شکل (۱7) تغییرات زمانی عمق آب‌شستگی برای مدل 
پایه‌های گروه ۴ در بازة زمانی ۱ ساعته نمایش داده شذه 
است. نتایج حاصل‌شده نیز در جدول (1) ارائه گردیده 
است. 

نتایج مدل عددی نشان می‌دهد پایة سکویی ٩‏ 
سانتی‌متری با کاهش ۲0-4690 بیشترین تأثیر را در کاهش 
فرسایتن بستتر داشعه است. این سکوی در بالادست. پابه 


نشریه مهندسی عمران فردوسی 


۳ 


مانع از برخورد جت جریان به بستر بالادست پایه می- 
گردد. مشهود است که با افزایش ارتفاع سکو از نقش 
محافظتی آن کاسته می‌گردد طوری که کاهش ۱۲ 
درصدی عمق آب‌شستگی مربوط به مدل سکویی ۸ 


جدول (1): درصد کاهش عمق آب‌شستگی در مدل پایه‌های گروه 


3 
ی ز و نم هم اس | ی 
0 ۰ ۰ ۳۹ ۰۸/۵ تکیت 
0591 ۶1۰ 5 ۳/۳۲ ۰/5۰۵ ۳۰ 
02 3 1 ۱۹/۹ ۰۳۹ 3 
0103 ۰ ۹ ۳۳ ۰/۷۵ ۳۱ 
54 ۰ ۹ ۳۹/۷ ۰,۲ ۱۲ 


در استفاده از ابزار سکو, تنها قسمتی از جریان رو 
به پایین در بالای سکو حذف گردیده و جریان رو به 
پایین در مقطع سکو دوباره به بستر بالادست برخورد 
کرده و باعث تشدید آب‌شستگی می‌گردد. با افزایش 
ارتفاع سکو بر شدات جت جریان رو به پایین در بالادست 
آن افزوده می‌شود. همچنین با بالا رفتن ارتفاع سکو 
درصد مساحت انسداد افزایش يافته و بدین ترتیب از 


عملکرد اصلی سکو کاسته می‌شود. در مدل‌های 083 و 


2000 2500 3000 3500 4000 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


4 تمایل به برقراری تعادل در شرایط فرسایش دیده 
نمی‌شود؛ اما در دو مدل سکویی ۵ و ٩‏ سانتی‌متری 
مجانب افقی در زمان ۳۲۰۰ ثانیه قابل مشاهده است. 
فرسایش بستر پیرامون مدل 0382 در شکل (۱۷) نشان داده 


یله اننتت: 


مدل‌های گروه ) 
شکل (۱۸) تغییرات زمانی عمق آب‌شستگی برای مدل- 
های گروه » را نمایش می‌دهد. بر حسب نتایج. ارتفاع 
مر سکو ۰ سانتی‌متر در نظر گرفته شده است. ترکیب 
مقطع هندسی با سکو به‌منظور بررسی میزان تأثیر هرکدام 
از پارامترهای آب‌شستگی» جدایش جریان و جریان رو 
به پایین؛ بوده است. جدول (۷) مقایسة کمّی بیشینة عمق 


آب‌شستگی را نشان می‌دهد. 


جدول (۷): درصد کاهش عمق شتا در مدل پایه‌های گروه 
۵ 
| رن | هم | چم اسص | 
0 ۰ ۰ ۳۹ ۰/۱/۹۵ سس 
9 ۶۰ 5 ۱9/۵ ۰/۱۳6۸ 3 
۹:1 ۰ ۹ ۳۹ ۰:۳۳ ۱۰ 
1۳۳ ۰ 1 ۱۸/۳ ۰/۳۵ ۳۷ 
11۳ ۶۰ 1 ۱۱ ۰۱۸۳ نژ 


(مزظ م1م۳ن6)) بعزظ 0عاعما۲0مرصلا 


م۳ تهلناوصم»۲6 11090-۲0560 ۰ 


۵۲ [هعا60[[ظ ---- 
6۲ مهلنم‌تاطمیا ‏ 


تگی (0۲0) 


0 500 1000 1500 


شکل (۱۵): تغییرات زمانی بیشینهٌ عمق آب‌شستگی در مدل‌های گروه ۸ 


نشريه مهندسی عمران فردوسی 


سال سی و سوم شمارة دی ۱۳۹۹ 


نازیلا کاردا - حبیب حکیم زاده- پوسف حسن زاده 


4000 3500 32000 2500 2000 
زمان (ثانیه) 


رهز 6۳۵6) معط 0ع۲ع۲۵1۵مرصلا 
۷۵061 052 
۷۵061 31 
006 83:) --. 
۷۵061 254 


عمق ی 


تکی () 


0 500 1000 1500 


شکل (۱1): تغییرات زمانی بیشینةٌ عمق آب‌شستگی برای مدل پایه‌های گروه 7 


(پ 


شک (۱۷):فرسایش ستر برامون مدل پاية 69 (الف) بردارهای سرعت جریان. کر سرآمون بای سکوتی لب) فرسایش یک‌ساه یستر 
پیرامون پايةُ سکویی (پ) نمایی سه‌بعدی از آب‌شستگی اطراف مدل 


با افزایش قطر پایه در مقطعی از آن عرض پایه در 
محل بستر افزايش می‌یابد که درنهایت موجب ازدیاد 
عرض انسداد می‌شود. ایجاد سکو و ازدیاد عرض انسداد 
در مدل‌های 18 و 8 موجب افزايش عمق آب‌شستگی در 
مدل‌های 1 و (شته است: دو مثل با ند تایر سکر 
تنها کاهش ۱ میلی‌متری عمق فرسایش بوده است. اگرچه 
افزایش قطر پایه در نزدیکی بستر مانع برخورد جت 
جریان به بستر می‌گردد؛ لیکن افزایش قطر باعث تشدید 
جدایش جریان در کناره‌های پایه شده و بر شدّت 


کی ی اف انلن 
مدل‌های گروه ۱ 
هدف استفاده از شکاف میانی» منحرف نمودن جریان رو 


به پایین در مقطع بالادست پایه و کاهش عرض انسداد 


سال سی و سوم. شمارهُ دو, ۱۳۹۹ 


است. عرض. ارتفاع و محل قرارگیری شکاف و طولی از 
شکاف که درون بستر رودخانه قرار می‌گیرد عوامل مهم 
در شکاف پایهٌ پلها هستند. در این تحقیق عرض شکاف 
مطابق با پیشنهادهای پژوهشگران مختلف. برابر با ۲۵ 
درصد قطر پایه و ارتفاع آن یک برابر قطر پایه انتخاب 
شده است که عرض شکاف برابر ۲ سانتی‌متر و ارتفاع 
ان ۵,سات‌ستر بلاق سس است, شکاف‌ها دی زیر راز 
بستر نیز امتداد یافته‌اند. با توجه به این‌که ارتفاع موثر 
سکو در مدل‌های گروه 2 1 سانتی‌متر تعیین شده است. 
ارتفاع شکاف ۱ سانتی‌متر کوچک‌تر از ارتفاع سکو در 
نظر گرفته می‌شود. مقایسه شماتیک تغییرات بيشینة عمق 
آب‌شستگی برای مدل‌های گروه ‏ در شکل (۱۹) ارائه 


میله: ابیت 


نشریه مهندسی عمران فردوسی 


2000 250 3000 3500 000 


زمان (ثانیه) 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


رد۳۱6 6۱۲۵۱6)) ۲عزظ ۲۵1۵0۲60 مرصل] 
6 3 
[۷۵06 ۳ .-.- 
۷0061 5 -- 


ات 


0 500 1000 1500 


شکل (۱۸): تغییرات زمانی بيشینة عمق آب‌شستگی برای مدل پایه‌های گروه » 


4000 


3500 


2500 3000 


(6۲ز۳ 6۳۵16)) عزظ 0ع6۲ع۲0]۲صط] 
ادا 
۷۵6 1185 ---- 
[ع۷0۵ 385 - - 


نتگی (701) 


1500 


شکل (۱۹): تغییرات بيشينة عمق آب‌شستگی برای مدل‌های گروه 1 


در لحظات آغازین شبیه‌سازی» عمق آب‌شستگی 
مدل‌های 188 ۴8 و 18 دارای تغییرات قابل توجهی 
است لیکن پس از گذشت زمان ۲۵۰۰ ثانیه این تغییرات 
کاهش می‌یابد تا اینکه درنهایت به حالت تعادل برسد. در 
مدل پایه‌های شکاف‌دار تغییرات اولیه در عمق 
آب‌شستگی بسیار بیشتر از تغییرات عمق آب‌شستگی در 
مدل پایه‌های گروه ۵ 3 و 6 است که نشان می‌دهد 
مدّت زمان لازم برای رسیدن به شرایط تعادل در مدل‌های 
شکاف‌دار بیش از ساير مدل‌ها بوده و زمان بیشتری برای 
شبیه‌سازی نیاز است. 

اعمال بازشدگی در مدل 6 تغییر قابل توجهی در 
عبقی آن‌شسشکی اد نگرده ات نات زیت شکافت 
در کاهش فرسایش بستر تنها در لحظات اولیه شبیه‌سازی 
دیده می‌شود لیکن با گذشت زمان (۲۷۵۰ تانیه 4 رسویات 


فرسایش يافته از کناره‌های پایه (ناشی از جدايش جریان) 


نشريه مهندسی عمران فردوسی 


به سمت پاییندست کشیده شده و در پشت پایه تجمع 
می‌یابند. تجمع رسوبات موجب گرفتگی شکاف شده و 
باعث می‌شود مدل ث) در بلندمدات رفتاری همانند مدل 
داشته باشد. پیشیتة عمق آب‌شستکین برای فدل خ) 
حدود ۲۷ میلی‌متر است که اندکی کمتر از بيشینة عمق 
آب‌شستگی مدل 6 است» اسبداد شکاف در مدل‌های 
8 28 و 18 مشاهده نمی‌شود. در این مدل‌هاء بزرگ 
بودن نسبت طول به عرض پایه موجب می‌شود تا با 
کاهش سرعت جریان در کناره‌های پایه. رسوبات 
برخاسته در کناره‌های مدل‌ها ته‌نشین شده و به سمت 
پایین‌دست جابه‌جا نشوند؛ بنابراین شکاف ابزار مناسبی 
برای کاهش فرسایش بستر در پایه‌هایی است که نسبت 
طول به عرض آن‌ها بزرگ‌تر از ۱ باشد. عمق فرسایش 
در مدل 185 ۱2/۵ میلی‌متر به دست آمده است که برابر 


۰ درصد کاهش نسبت به مدل مبنا و ۳۱ درصد کاهش 


سال سی و سوم شمارة دی ۱۳۹۹ 


نا زیلا کاردان- حبیب حکیم زاده- پوسف حسن زاده 


نسبت به مدل 3 است. 

به عبارتی ترکیب مقطع مستطیلی گرد گوشه با سکو 
و شکاف باعث کاهش عمق آب‌شستگی به میزان ۱۳/۵ 
میلی‌متر شده است. در مدل 85 کاهش بيشينة عمق 
آب‌شستگی نسبت به مدل مبنا و مدل * به ترتیب برابر 
۸ درصد و 1 درصد است. اعمال سکو و شکاف در 
مدل بآ تنها موجب کاهش ۱ میلی‌متر در بيشينة مقدار 
فرسایش بستر شده است. در مدل 8٩‏ نیز بيشینة عمق 
آب‌شستگی به‌اندازه 7۲ درصد نسبت به مدل مبنا و ۶ 
درصد نسبت به مدل بآ کاهش داشته است. این مقدار 
کاهش در مقایسه با تأثیر سکو و شکاف در مقطع 
مبتخطیلی گرد گنه ییاز انلدک, استه 

در مدل با مقطع هندسی بیضی و عدسی با ایجاد 
شکاف. شکل مقطع هندسی تغییر کرده و پایه دیگر 
آیرودینامیک نبوده و از انطباق آن با الگوی جریان کاسته 
می‌شود؛ بنابراین مقطع هندسی پایه تأثیر اندکی بر کاهش 
آب‌شستگی خواهد داشت. در مدل پایه .1 عمق فرسایش 
بستر در لحظات آغازین شبیه‌سازی بسیار کم است لیکن 
در مدل 1.3 از لحظات آغازین اجرای برنامه. تشدید 
فرسایش در اطراف لبه‌های شکاف دیده می‌شود. 
به‌عنوان‌مثال عمق فرسایش در دقیقه ۲۰ شبیه‌سازی برای 
مدل بآ برابر ۰/۵ میلی‌متر است درحالی‌که این مقدار در 
مدل 18 مساوی ۱/۵ میلی‌متر است؛ بنابراین تأثیر مثبت 
ایجاد شکاف در مقطع هندسی عدسی به‌مراتب کمتر از 
تأثیر آن در مدل با مقطع مستطیلی گرد گوشه است. 
شکل‌های (۲۰) و (۲۱) نمای دوبعدی و سه‌بعدی از 
آب‌شستگی اطراف پایه‌های مرکب 2۳858 و 185 را 


«> ۱ 
۱۱۱ 
۹۹۰۹ ۳۷۴۳۹۳۱۱ 


اف (ب) 


شکل (۲۰): فرسایش بستر پیرامون مدل 185 


سال سی و سوم. شمارهُ دو, ۱۳۹۹ 


(الف) (ب) 


در جدول (۸) درصد کاهش بیشینة عمق 
آنب‌شستکین برای, مدل بانه‌های. گروه ۳ نیت به هدن 


مبنا نشان داده شده است. 


جدول ۸ بيشینةٌ عمق آب‌شستگی مدل پایه‌های گروه ط 


خی | و | هم | هم | سصث | مق 
0 ۰ 3 ۳۹ ۰۱/۵ سس 
0155 ۰ 1 ۳۷ ۵ 1/۹ 
دص ۰ 1 ۱/۵ ۱:۱" 0۰ 
۳۳95 ۰ 1 ۱۵ ۰/۳6۵ 1۸ 
دص ۰ 1 ۱۱ 12 ن 


عملکرد مدل پایه‌های گروه ۸ در کنترل فرسایش 
بستر مطلوب‌تر از مدل‌های گروه 8 بوده است. یعنی 
استفاده از شکل هندسی آیرودینامیک در پایه‌های پل. 
نقش حفاظتی موثری در کاهش فرسایش بستر داشته 
است. در مدل‌های سکودار اگرچه ایجاد سکو در پایه 
موجب کنترل جریان‌های رو به پایین و گردابه‌های 
نعل‌اسبی می‌گردد لیکن افزایش قطر پایه در تراز بست 
موجب افزايش عرض انسداد و تشدید فرسایش بستر 
می‌شود. لذا کتترل جدایش جریان بسیار موثرتر از کنترل 
جریان رو به پایین در کاهش آب‌شستگی است. ترکیب 
سکو, شکل پایه و شکاف میانی نیز کاهش قابل توجه در 
عمق آب‌شستگی را موجب شده است لیکن در اين پایه‌ها 
احتمال گرفتگی و انسداد شکاف در اثر جابه‌جایی 
زسوزباث جرک دازد فر این پایه ابا مادک کاف: سم 
تایه کاتی بای کاس سای زا بایان 


نشریه مهندسی عمران فردوسی 


دارد؛ بنابراین استفاده از پایه‌هایی با مقاطع آیرودینامیک 


می‌تواند به‌عنوان پيشنهاد کارهای آتی مطرح گردد. 


نتیجه گیری 

در پژوهش حاضر شبیه‌سازی عددی پایه‌های پل ساده و 
مرکب با استفاده از نرم‌افزار ۳10380 انجام گرفته است. 
به همین منظور در ابتدا گزارشی از نحوه مدل‌سازی پایه- 
های پل در شرایط جریان زلال رودخانه‌ای ارائه شد که 
شامل مشخصات مدل فیزیکی و مدل عددی است. در 
ادامه, صت‌سنجی پارامترهای مورد استفاده مانند اندازه 
سلول‌های شبکه مدل‌های آشفتگی و عمق فرسایش 
مورد بررسی قرار گرفت. 

مدل پایه‌های مرکب با هدف بررسی دو پارامتر مهم 
در وقوع آب‌شستگی (جت جریان رو به پايین و جدایش 
جریان) مورد مطالعه قرار گرفته و درصد کاهش عمق 
آب‌شستگی آن‌ها نسبت به مدل پایه استوانه‌ای به دست 
آمد. نتایج نشان داد در اغلب موارد استفاده از یک ابزار 
کنترل آب‌شستگی بهتر از ترکیب ابزارهای مختلف عمل 
می‌نماید. به‌کار گیری شکل عدسی به‌عنوان مقطع هندسی 
پایه» توانست عمق آب‌شستگی را تا ۵۸ درصد کاهش 
دهد. مدل ترکیبی مقطع مستطیلی گردگوشه. سکو و 
شکاف نیز توانست فرسایش بستر را تا ۵۰ درصد کاهش 
دهد لیکن استفاده از مدل پایه‌های شکاف‌دار برای 


ترا 


شبیه‌سازی سه‌بعدی پدیده رهایی متناوب گردابه‌ها و انتقال رسوب بستر... 


رودخانه‌هایی توصیّه می‌شود که در آنها حجم جابه‌جایی 
رسوبات پایین باشد. 


سرعت متوسط 10 
فشار ۲ 
مساحت ۸ 
شتاب ناشی از نیروهای وزنی 0 
شتاب ناشی از نیروهای ویسکوز 1 
ی برشنی فیوازد 
دانسیته سیّال 0 
ثابت ون کارمن تا 
ضریب زبری 1 
کسر مربوط به قسمت سیّال َ 
سرعت شناوری ذرات رسوب منوا 
ضریب درگ ذرّه رسوب 1 
ضریب شناوری ۷ 
قطر متوسط ذرات رسوب ی 
تنش برشی بحرانی 3 
پارامتر شیلدز بحرانی 9 
سرعت بلند شدن رسوبات زرا 
بردار نرمال 1 
زاوية اصطکاک داخلی رسوبات 3 
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نشريه مهندسی عمران فردوسی 
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فقم رنه ناس مرن فزگارسین سال سی و سوم شماره دو ۱۳۹۹ 


